Theory of charge transport in diffusive normal metal / conventional superconductor 
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Tunneling conductance in diffusive normal metal / insulator / s-wave superconductor (DN/I/S) 
junctions is calculated for various situations by changing the magnitudes of the resistance and 
Thouless energy in DN and the transparency of the insulating barrier. The generalized boundary 
condition introduced by Nazarov [Superlattices and Microstructures 25, 1221 (1999)] is applied, 
where the ballistic theory by Blonder, Tinkham and Klapwijk (BTK) and the diffusive theory 
by Volkov, Zaitsev and Klapwijk based on the boundary condition of Kupriyanov and Lukichev 
(KL) are naturally reproduced. It is shown that the proximity effect can enhance (reduce) the 
tunneling conductance for junctions with a low (high) transparency. A wide variety of dependencies 
of tunneling conductance on voltage bias is demonstrated including a )7-shaped gap like structure, 
a zero bias conductance peak (ZBCP) and a zero bias conductance dip (ZBCD). The temperature 
dependence of tunneling conductance is also calculated and the conditions for the reentrance effect 
are studied. 



I. INTRODUCTION 

The electron coherence in mesoscopic superconduct- 
ing systems is one of the important topics of solid state 
physics. The low energy transport in these systems 
is essentially influenced by the Andreev reflection^, a 
unique process specific for normal metal/superconductor 
interfaces. The phase coherence between incoming elec- 
trons and Andreev reflected holes persists in the diffu- 
sive normal metal at a mesoscopic length scale and re- 
sults in strong interference effects on the probability of 
Andreev reflection^. These effects become prominent at 
sufficiently low temperatures where the thermal broad- 
ening is negligible. One of the remarkable experimen- 
tal manifestations is the zero bias conductance peak 
(ZBCP)3iWi7i8i9ii°i"ii2ii3^ A calculation of tunneling 
conductance in a normal metal (N) / superconductor (S) 
junction is an interesting theoretical problem since quan- 
tum interference effects due to Andreev reflection are ex- 
pected. 

For a clean NS contact in the presence of the inter- 
face potential barrier the conductance was calculated by 
Blonder, Tinkham and Klapwijkii (BTK) in terms of 
the corresponding transmission coefficients on the basis 
of the solution of the Bogoljubov - de Gennes equations. 
From the general set of boundary conditions connect- 
ing the quasiclassical Green's functions on both sides 
of the interface for arbitrary transmission probabilities, 
Zaitsev^^ derived the expression for the conductance sim- 
ilar to that by BTK. The BTK methodic is confined to 
ballistic systems. The generalization of this method to 
systems with impurities has been performed by several 
authors (see the reviewi^). In a number of papers the 
transmission coefficients were directly calculated by nu- 
merical methods^^'^^. However, it is difficult to apply 
such methods to most of relevant experimental situa- 
tions. Another approach, the so-called random matrix 



theory, was employed by Beenakker et al., where the 
total transmission coefficients are expressed in terms of 
those through the normal part of the system and the 
normal/superconductor interface separateljsiSiiS. Within 
this theory, the ZBCP observed in experiments is under- 
stood as a resonance phenomenon related to reflectionless 
tunneling^*^ . The scattering matrix approach was later 
generalized to finite voltage and temperature^ 

On the other hand, a quasiclassical Green's function 
calculation based on nonequilibrium superconductivity 
theories^^ is much more powerful and convenient for the 
actual calculations^'^. In this approach, the impurity 
scattering is included in the self-consistent Born approx- 
imation and the weak localization effects are neglected. 
In the theory of tunneling conductance developed by 
Volkov, Zaitsev and Klapwijk (VZK) by solving the Us- 
adel equationsSi, the origin of the ZBCP observed in sev- 
eral experiments was clarified to be due to the enhance- 
ment of the pair amplitude in the diffusive normal metal 
by the proximity effect. VZK applied the Kupriyanov 
and Lukichev (KL) boundary condition for the Keldysh- 
Nambu Green's function^^. The KL boundary condition 
is valid for the atomically sharp interface barrier dividing 
two diffusive metals. As shown by Lambert et alm^, this 
condition is exact in two limits of either high or low bar- 
rier transparency, with small corrections in the interme- 
diate transparency range. By applying the VZK theory, 
several authors studied the charge transport in various 



junctions 



27.28.29.30.31.32.33.34 



by solving the Usadel equa- 



tions, (see review by Belzig et alMi). 

The generalization of the KL boundary conditions for 
an arbitrary connector between diffusive metals was pro- 
vided by Nazarov within the so-called "generalized cir- 
cuit theory"J^. In this theory, the mesoscopic system 
is presented as a network of nodes and connectors. A 
connector is characterized by a set of transmission coef- 
ficients and can present anything from a ballistic point 
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contact to a tunnel junction. A conservation law of ma- 
trix current holds in each node. The method to derive the 
relation between matrix current and Green's functions 
puts the results of Ref.-'^^ to the framework of Landauer- 
Biittiker scattering formalism. The boundary condition 
for Keldysh-Nambu Green's function was derived in2& for 
an arbitrary connector including various situations from 
ballistic point contact to diffusive contact. Actually, this 
boundary condition is very general since the BTK theory 
is reproduced in the ballistic limit while in the diffusive 
limit with a low transmissivity of the interface, the KL 
boundary condition is reproduced. 

Although a number of papers were published on charge 
transport in mesoscopic NS junctions, as far as we know, 
almost all of them are either based on the KL bound- 
ary conditions or on the BTK model. However in the 
actual junctions, transparency of the junction is not nec- 
essarily small and impurity scattering in the DN is im- 
portant. Therefore, an interesting and important theo- 
retical problem is the calculation of the tunneling con- 
ductance in normal metal / conventional superconduc- 
tor junctions using the boundary condition from Ref.'^^ 
since both the ballistic (BTK theory) and diffusive (VZK 
theory) cases can be covered simultaneously. In the 
present paper, we study the tunneling conductance in 
diffusive normal metal / insulator / conventional su- 
perconductor (DN/I/S) junctions for various parameters 
such as the height of the insulating barrier at the in- 
terface, resistance Rd in DN and the Thouless energy 
Exh in DN. We concentrate on the normalized tunnel- 
ing conductance of the junctions arieV) as a function 
of the bias voltage V. The conductance arieV) is given 
by arieV) — <7s{eV)/aN{eV) where as[N){eV) is the 
tunneling conductance in the superconducting (normal ) 
state at a bias voltage V. 

In the present paper the following points are clarified: 

1. When the transparency of the junction is sufficiently 
low, aT{eV) for | eT^ |< Aq is enhanced with the in- 
crease of Rd due to the enhancement of the proximity 
effect. The ZBCP becomes prominent for Eth « Ao 
and Rd/Rb < 1. In such a case, with a further increase of 
Rd/Rb the ZBCP changes into a zero bias conductance 
dip (ZBCD). In the low transparent limit, the line shapes 
of aT{eV) are qualitatively the same as those obtained 
by VZK theorjiSiSa. 

2. When the transparency of the junction is almost unity, 
axieV) always exhibits a ZBCD except for the special 
case of Rd = 0, i.e. the BTK limit. 

3. The measure of the proximity effect, 9, is mainly 
determined by Rd/Rb and Eth, where Rb is the resis- 
tance from the insulating barrier. The proximity effect 
enhances (reduces) the magnitude of axieV) for junc- 
tions with low (high) transparency. 

4. Even for junctions between conventional s-wave super- 
conductors, we can expect a wide variety of line shapes 
of the tunneling conductance, a ZBCP, ZBCD, [/-shaped 
structure, and a rounded bottom structure. 

5. We have clarified the parameter space where a ZBCP 



should be expected. Typically, small Thouless energy 
Erh is required for a ZBCP. If the magnitude of Exh is 
increasing up to Aq, a ZBCP is only expected for junc- 
tions with low transmissivity, Rd/Rb << 1- 

The organization of this paper is as follows. In section 
2, we will provide the detailed derivation of the expression 
for the normalized tunneling conductance. In section 3, 
the results of calculations of arieV) and 6 are presented 
for various types of junctions. In section 4, the summary 
of the obtained results is given. 



II. FORMULATION 

In this section we introduce the model and the for- 
malism. We consider a junction consisting of normal 
and superconducting reservoirs connected by a quasi-one- 
dimensional diffusive conductor (DN) with a length L 
much larger than the mean free path. The interface be- 
tween the DN conductor and the S electrode has a resis- 
tance Rd while the DN/N interface has zero resistance. 
The positions of the DN/N interface of the DN/S inter- 
face are denoted as a; = and x — L, respectively. Ac- 
cording to the circuit theory, the interface between DN 
and S is subdivided into two isotropization zones in DN 
and S, two ballistic zones and a scattering zone. The sizes 
of the ballistic and scattering zones in the current flow 
direction are much shorter than the coherence length. 
Although the generalized boundary condition of Refi^ 
is valid for arbitrary interfaces, here scattering zone is 
modelled as an infinitely narrow insulating barrier de- 
scribed by the delta function U{x) = Hd{x — L). The 
resulting transparency of the junctions Tm is given by 
r,„ = 4cos2 (/)/(4cos2 (/) + Z^), where Z = iH/iTivp) is a 
dimcnsionless constant, (j) is the injection angle measured 
from the interface normal to the junction and vp is Fermi 
velocity. Variation of the barrier shape will not change 
our results in the considered case of isotropic supercon- 
ductivity in the S electrode. 

We apply the quasiclassical Keldysh formalism in the 
following calculation of the tunneling conductance. The 
4x4 Green's functions in DN and S are denoted by 
Gi(x) and G2{x) which are expressed in matrix form as 



Gi{x) = 
G2{x) = 



Ri{x) Kiix) 

Ai{x) 

i?2(x) k2{x) 

A2{x) 



(1) 

(2) 



where the Keldysh component Ki^2{x) is given by 
i^i(2)(a;) = ^i(2)(a;)/i(2)(a;) - /i(2) (2^)^1(2) (a^) with re- 
tarded component -Ri.2(a;), advanced component ^1^2(2^) 
using distribution function fi{2){x)- In the above, R2{x) 
is expressed by 



3 



Ao/v/A§"-e2 



with g = e/i/e^ - and / - y 
where e denotes the quasiparticle energy measured from 
the Fermi energy, A2{x) = — i?2(a;) and f2{x) = 
tanh[e/(2kBT)] in thermal equihbrium with temperature 
T. We put the electrical potential zero in the S-electrode. 
In this case the spatial dependence of Gi {x) in DN is de- 
termined by the static Usadel equationSi, 



D^[G,{x)^^^]+i[H,G,{x)]=Q, (3) 

with the diffusion constant D in DN, where H is given 
by 



H = 



i?o 

^0 y ' 



with Hq — er^. 

The boundary condition for Gi{x) at the DN/S inter- 
face is given by Nazarov's generalized boundary condi- 
tion. 



(4) 



B 



2T„[Gi(L_),G2(£+)] 
4 + T™([Gi(L_),G2(£+)]+-2)' 



The average over the various angles of injected particles 
at the interface is defined as 



< / dcf) cos (t)B{(t>) I i d(t>T{4>) cos (f) 

J--IT/2 J-tt/2 



with B{4>) = B and T{4>) = T„. The resistance of the 
interface Rh is given by 

9 

Rb = Rq 



rf0r(0)cos0' 

Here Rq is Sharvin resistance, which in three-dimensional 
case is given by R^^ — e'^kpSc/ {4:n^fi), where kp is the 
Fermi wave-vector and Sc is the constriction area. Note 
that the area Sc is in general not equal to the crossec- 
tion area Sd of the normal conductor, therefore Sc/ Sd is 
independent parameter of our theory. 

For Tm ^ in Eq. (4), the quantity B can be expressed 

as 

and we can reproduce the KL boundary condition. On 
the other hand, at a; = Gi{Q) coincides with that in the 
normal state. 

The electric current is expressed using Gi{x) as 



4ei?, 



d Jo 



'deTr[r3(Gi(x)^^)^], (5) 



where (Gi(x) "^q!^'' denotes the Keldysh component 



of (Gi(a;)^^^j^). In the actual calculation it is conve- 
nient to use the standard ^-parameterization when func- 
tion Ri{x) is expressed as 



Ri{x) = fs cosO{x) + T2 sm9{x). 



(6) 



The parameter 9{x) is a measure of the proximity effect 
in DN. ^ ^ 

Functions Ai{x) and Ki{x) are expressed as Ai{x) = 
-Rl{x) and Ki{x) = Ri{x)fi{x) - fi{x)Ai{x) with the 
distribution function fi{x) which is given by fi{x) = 
fi{x) +T3ft{x). In the above, ft{x) is the relevant distri- 
bution function which determines the conductance of the 
junction we are now concentrating on. From the retarded 
or advanced component of the Usadel equation, the spa- 
tial dependence of 9{x) is determined by the following 
equation 
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D—e{x) + 2iesin[e{x)] = 0, 
while for the Keldysh component we obtain 

^^[^cosh2^.™..(^)]=0. 



(7) 



(8) 



At a; = 0, since DN is attached to the normal electrode, 
6'(0)=0 and /t(0) = /to is satisfied with 

/to = \{ianh[{e+eV)/{2kBT)]-t&nh[{e~eV)/{2kBT)]}. 

Next we focus on the boundary condition at the DN/S 
interface. Taking the retarded part of Eq. Q), we obtain 



L de{x) 
Rd dx 



<F> 



|x=L_ = 



Rh 



(9) 



F = 



(2- 



2{f cos OL-g sin eL)T,n 

Tm) +T^[g COS 9l + fsindLY 



with Ol = 0{L-). 

On the other hand, from the Keldysh part of Eq. 10} 
we obtain 

'R^ ~dx ^imag{x) \x=L^= 



Rb 



with 



(10) 



T,2„Ai-f2T„,(2-T„0A2 
''° 2 I (2 - r„, ) + T„ [g COS Bl + J sin Ol] 



12 ' 



Ai - (1+ I cosBl P + I sinOL Hd 9? + \f? +1) 



4Imag[/g*]Imag[cos sin 6*jJ\ 



(11) 
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A2 = ReallgCcosfiiL + cos6'l) + f(sin6iL + sin^iL)}, (12) 

where Oimag{x) denotes the imaginary part of 6{x). For 
Tm « 1, Ibo is reduced to 



'60 



Real[g(cos + cos 91) + f(sin Ol + sin 61)] , 



(13) 

which is the expression used in the VZK theory. After a 
simple manipulation, we can obtain ft{L-) 



ML-) = 



COSh-^ Oimag {x) 



Since the electric current lei can be expressed via 0i in 
the following form 

L d ft 

lei = 5- / (^) \x=L. cosh^[Imag(6'i)]de, 

eRd Jo ox 

we obtain the following final result for the current 



Iel = - 



de- 



fi 



to 



Rb 4_ Rd 

<ho> L Jo 



dx 



(14) 



cosh^ 9imag{x) 



Then the total resistance R at zero temperature is given 
by 



R = 



Rb 



Rd 



< ho > 



dx 



cosh 9i. 



,{x) 



(15) 



and the tunneling conductance in the superconducting 
state (Ts{(iV) is given by as(ey) = 1/R. 

It should be mentioned that for Rd = 0, Ol becomes 
zero due to the absence of the proximity effect. Then Ibo 
is given as follows 



-^60 = 



(1+ I 5 I' + I / \')T^ + 2T„(2 - T„)Real(5) 



(2 - Tm) + Tmg 



-(Tm - 1) I r n 



1 - (1 - T„)r2 



(16) 



with r 



and the resulting as is given by 



ag{eV) = —l cos (j)d(j), 

and reproduces that by BTK theory. 

It should be remarked that in the present circuit the- 
ory, Rd/ Rb can be varied independently of T^, i.e. inde- 
pendently of Z, since one can change the magnitude of 
the constriction area 5c independently. In other words, 
Rd/Rb is no more proportional to Tav{L/l), where Tav 
is the averaged transmissivity of the barrier and I is the 



mean free path in the diffusive region, respectively. Based 
on this fact, we can choose Rd/Rb and Z as independent 
parameters. 

In the following section, we will discuss the normalized 
tunneling conductance arieV) = as{eV)/aN{eV) where 
(Tjv(eV^) is the tunneling conductance in the normal state 
given by crjv(eF) = cfn = ^/{Rd + Rb), respectively. 



III. RESULTS 



Tunneling conductance vs voltage: zero-bias 
anomalies 



In this section, we focus on the line shape of the tunnel- 
ing conductance. Let us first choose the relatively strong 
barrier Z = 10 (Figs. 1 and 2) for various Rd/Rb- For 
Exh = Ao, the magnitude of arieV) for | eV |< Ag in- 
creases with the increase of Rd/Rb- First, the line shape 
of the tunneling conductance remains to be U shaped and 
only the height of the bottom value is enhanced (curve h 
). Then, with a further increase of Rd/Rb, a rounded 
bottom structure (curve c and d) appears and finally 
it changes into a nearly fiat line shape (curve e). For 




eV/Ao 



FIG. 1: Normalized tunneling conductance for Z=10 and 
Eth/Ao = 1. a: Rd/Rb = 0, b: Ra/Rb = 0.1, c: Ra/Rt = 1, 
d: Rd/Rb = 2 and e: Rd/Rb = 10. 

Exh = O.OlAo (Fig. 2), the magnitude of arieV) has 
a ZBCP once the magnitude of Rd/Rb deviates slightly 
from 0. The order of magnitude of the ZBCP width is 
given by Eth- When the magnitude of Rd/Rb exceeds 
unity, the ZBCP splits into two (curve d) and finally 
axi&V) acquires a zero bias conductance dip (ZBCD) 
(curve e). In the limit of an extremely strong barrier 
with Z ^ QG, the results of the VZK theory are repro- 
duced. 

On the other hand, in the fully transparent case with 
Z = the line shape of the tunneling conductance be- 
comes quite different. For Eth = Aq (Fig. 3), the mag- 
nitude of axieV) decreases with the increase of the mag- 
nitude of Rd/Rb- The bottom parts of all curves are 
rounded and arieV) always exceeds unity. On the other 
hand, for Exh = O.OIAq, (TT{eV) has a ZBCD even for a 
small magnitude of Rd/Rb except for the special case of 
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Rd/Rb — where the BTK theory is vahd. This feature 
is quite different from that shown in Fig. 2. In the case 
of an intermediate barrier strength, Z = 1, the shape of 
arieV) becomes rather complex. For Eth = ^o, (TT{eV) 
has a shallow gap structure similar to the case of the BTK 
theory (curve a ). With the increase of Rd/Rb. the co- 
herent pc^ak structure at eV = ±Ao is smeared out and 
the voltage dependence becomes very weak as shown by 
curve e. For Eth — O.OIAq, the ZBCP appears for a 
small magnitude of Rd/Rb (see curve 6). With increas- 
ing magnitude of Rd/Rb, the ZBCP changes into a ZBCD 
(curves c d e). As compared to the Z = 10 case (see Fig. 
2), the ZBCP is much more easier to change into a ZBCD 
by increasing the magnitude of Rd/Rb- 

It is interesting to study how various parameters influ- 
ence the proximity effect. The measure of the proximity 
effect at the S/N interface 6l is plotted for Z = and 
Z = 10 with corresponding parameters in Figs. 1 to 4. 
For Rd/Rb = 0, 9l = is satisfied for any Eth and Z. 
Besides this fact, at e = 0, 9l always becomes a real 
number. First, we study the case of Eth/^o = 1 (Fig. 
7) where the same values of Rd/Rb are chosen as in Figs. 
1 and 3. The real part of 9l is enhanced with an increase 
in Rd/Rb and is almost constant as function of e. At 
the same time, the imaginary part of 6l is an increasing 
function of e for all cases. There is no clear qualitative 
difference between the energy dependencies of Real(^^i) 
and Imag(6'L) at Z — and Z = 10. Next, we discuss 
the line shapes of 9l for Eth/^o = 0.01. Real(^i,) has 
a peak at zero voltage and decreases with the increase 
of e. Imag(0L) increases sharply from and has a peak 
at about e --^ Eth, except for a sufficiently large value of 
Rd- Also in this case, there is no qualitative difference 
between the line shapes of Real(Imag)(6i,) for the Z = 
case and that for Z = 10. 

Although the magnitude of 9l, the measure of 
proximity effect, is enhanced with increasing Rd/Rb, its 
influence on arieV) is different for low and high trans- 
parent junctions. In the low transparent junctions, the 
increase in the magnitude oi 9l by Rd/Rb can enhance 
the conductance axieV) for eV ^ and produce a 




1 

eV/A„ 



FIG. 3: Normalized tunneling conductance for Z = and 
Eth/ Aq = 1. a: Ra/Rt = 0, b: Rd/Rb = 0.1, c: Ra/Rb = 1 
d: Rd/Rb = 2 and e: Rd/Rb = 10. 




FIG. 4: Normalized tunneling conductance for Z = and 
Eth/Aq = 0.01. a: Rd/Rb = 0, b: Rd/Rb = 0.1, c: Ra/Rb = 
1 d: Rd/Rb = 2 and e: Rd/Rb = 10. 



ZBCP, whereas in high transparent junctions the en- 
hancement of 9l suppresses the magnitude of (7T{eV). 

Finally, we focus on the condition where the ZBCP 
appears. We change Z and Rd/Rb for fixed Eth- An 
upper critical value of Rd — Rbu exists where the ZBCP 
vanishes for Rd > Rbu- For Eth = O.OIAq, Rbu increases 
with Z and converges at nearly 1.4. For Exh = O.SAq, 
the lower critical value of Rd = Rbi also appears where 
the ZBCP vanishes for Rd < Ru- The ZBCP is expected 




eV/A„ 




eV/Ao 



FIG. 2: Normalized tunneling conductance for Z—10 and 
Eth/ Ao = 0.01. a: Rd/Rb = 0, b: Rd/Rb = 0.1, c: Rd/Rb = 
1, d: Rd/Rb = 2 e: Rd/Rb = 10. 



FIG. 5: Normalized tunneling conductance for Z = 1 and 
Eth/ Ao = 1. a: Ra/Rb = 0, b: Rd/Rb = 0.1, c: Ra/Rb = 1 
d: Rd/Rb = 2 and e: Rd/Rb = 10. 
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FIG. 8: Real part of Qh (upper panels) and imaginary part 
of it (lower panels) is plotted as a function of e. Z—10 
(left panels) and Z—Q (right panels) with Eth/^o = 0.01. 
a: Rd/Rb = 0.1 b: Rd/Rb = 1, c: Rd/Rb = 2 and d: 
Rd/Rb = 10. 
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■° 1 
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l^'zBCP 


dzBCP 



10 



10 



20 



FIG. 7: Real part of 6l (upper panels) and imaginary part 
of it (lower panels) is plotted as a function of e. Z=10 
(left panels) and Z=0 (right panels) with Exh/^o = 1. 
a: Rd/Rb = 0.1 b: Rd/Rb = 1, c: Rd/Rb = 2 and d: 
Rd/Rb = 10. 



FIG. 9: The parameter space where ZBCP appears. 



tive of Ibo at e = is given by 



for Rhi < Rd < Rbu- The magnitude of Rbu is suppressed 
drastically as compared to that for Exh = O.OIAq. For 
Erh > Ao, a ZBCP region vanishes for < Z < 20. In 
order to understand the crossover from the ZBCP to the 
ZBCD in much more detail, it is interesting to calculate 
the second derivative of the total resistance i? as a func- 
tion ofe = eyate = 0. For simplicity, here we focus 
on the case of a sufficiently large Z. For e < Aq, for 
simplicity total resistance is written as 



R — Ri + R21 



(17) 



sin 6'L,r , d'^OL,r , ^ /d9L,i.2 



since = and O^^i = is satisfied at e = 0. 

The resulting second derivative of the total resistance 
at zero energy is given as 

d^R _ d^Ri a^i?2 
de^ de^ de^ ' 



d^Ri 



Ri 



Rb 



<IbO> 



R2 = ^ 



dx 



L Jq cosh^^j(a;)' 



< IbO >= fsmeL,rCOSh9L,i, 



' +COS6<L,r o +Sm6'l,,r(- 



de 



(18) 



with 6]^ = 0}^ r+i(^L,ij where 9^ ^ and 9^ i are the real and 
imaginary parts oi 9l, respectively. The second deriva- 



a2i?2 2i?d f ,99,{x)^2^^ 



de 



(19) 
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some cases. The order of 



where 9i{x) is the imaginary part of 9{x). The sign of 
becomes negative and it can induce > m 

and i^^f-y is propor- 
tional to the inverse of E^f^ . The sign of is crucially 
influenced by the relative magnitude of the second term 
at the right-hand side of Eq. H18|l . On the other hand, 

the sign of is always negative. 

For Exh << Aq and small magnitude of Rd/Rb, the 
resulting 9L,r is sufficiently small and the magnitude of 

becomes positive. When this positive contribution 

overcomes the negative contribution from we can 

expect a resistance minimum at zero energy, i.e., a ZBCP. 
However, with increasing Rd/Rb, the magnitude of 
decreases due to the enhancement of the third term in the 
right-hand side of Eq. 1)18(1 . while increases. Then, a 
critical value Rd = Rbu appears, above which the ZBCP 
changes into a ZBCD with an increase of Rd/Rb- This is 
the mechanism of the crossover from a ZBCP to a ZBCD. 

When the magnitude of Eth is enhanced, the mag- 
nitudes of the second and third terms in Eq. H18|) are 
reduced and the first term can not be neglected. The 
resulting magnitude of is suppressed and the value 
of Rbu is reduced. This is the origin of the difference in 
Rbu in Figs. 9(a) and (b). 

It is also interesting to present similar plots as Fig. 9 
using an averaged transparency of the junction Tav 



(20) 



The results are shown in Fig. 10. For Exh = O.OlAo, 
the magnitude of Rbu decreases monotonically with in- 
creasing Tav and it vanishes about Tav ~ 0.8. For 
Eth = O.8A0, the magnitude of Rbi increases for an in- 
creasing Tav while that for Rbu decreases. The ZBCP 
region is restricted to small values of Rd/Rb- 



the measure of the proximity effect 9l is enhanced for 
I e |< Eth and the zero bias conductance is proportional 
to Tav- However, the magnitude of 6*^ at finite energy 
in the range Exh <| e |< Aq is drastically suppressed as 
shown in Fig. 8, due to the existence of the proximity 
induced minigap^i in the normal diffusive part DN of the 
order of Eth- As a result, in this regime the conductance 
channel that provides a contribution proportional to Tav 
is not available. Thus only the low voltage conductance is 
enhanced. On the other hand, for large Ej^h with E^h ~ 
Ao, the measure of proximity effect 9l is insensitive to 
energy for | e |< Aq (see Fig. 7), then the resulting 
(Tt (eV) is always enhanced for | eV^ | < Aq and the degree 
of the prominent enhancement of (tt(0) is weakened. 



B. Tunneling conductance vs temperature: 
reentrance effect 



Finally, we look at temperature dependence of conduc- 
tance for various Z and Rd/Rb and focus on the relevance 
to the corresponding results in VZK theory based on the 
KL boundary condition. We calculate tunneling conduc- 
tance at non zero temperature following 



c7sieV,T) = dhi/dV. 



(21) 



Then we define deviation of tunneling conductance from 
that at zero temperature given by 



6as^cTsieV,T)-asieV,0). 



(22) 



In the following, we will plot Sas/crN as a function of 
temperature T. 



Z=10 



DC 
DC 1 



ETh=0.01Ao 


ETh=0-8Ao 




(b) 








ZBCP^ 



0.5 



0.1 



FIG. 10: Similar plots as in Fig. 9 using Tav- 

In any way, the preferred condition for the formation 
of a ZBCP is the combination of the low transparency 
of the junction and the smallness of the Eth/ ^0 ratio. 
This situation is understood as follows. It is well known 
from the BTK theory that the magnitude of the zero bias 
conductance is almost proportional to Tav for Rd — for 
Tav << 1- With the increase in the magnitude of Rd, 
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T/E„ 



FIG. 11: Sas/uN is plotted as a function of T. Z=10 and 
Eth = O.OlAo. a: Ra/Rb = 0.1 b: Rd/Rb = 1, c: Ra/R^ = 2 
and d: Rd/Rb = 10. 

For Rd/Rb = 0.1 (curve a) and Rd/Rb — 1 (curve b), 
due to the existence of ZBCP as shown in Fig. 2, 5as/<JN 
takes negative value and decreases with T. While for 
Rd/Rb — 2, 6as/cr]y first increases and decreases again 
(curve c). With the further increase of the magnitude 
of Rd/Rb, Sus/un becomes positive and it has maxi- 
mum at a certain temperature. This effect is known as 
" reentrance effect" . It was predicted theoretically within 
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the VZK theory in^Si^fli^ and observed experimentally^. 
According to the theorji22i^2i^, for Ri, ^ 0, the maxi- 
mum value of 5as/(TN is about 0.09 at temperature of 
the order of Thouless energy. Note that the reentrance 
of the metallic conductance occurs as a function of bias 
voltage as well, but here we concentrate on the tempera- 
ture dependence since it was studied in most detail within 
the VZK theory. 



Rd/Rb=iO 




FIG. 12: Sas/cFN is plotted as a function of T. Rd/Rb ~ 10 
and Erh = 0.01 Aq. a: VZK theory, b: Z = 10, c: Z = 1 and 
d: ^ = 0. 

In order to study the consequences of our theory for 
the reentrance effect, we have calculated the temperature 
dependence of Sas/a-N for various Z for fixed Rd/Rb- 
For the comparison with the standard VZK theory, we 
also plot Sas/a-N using the KL boundary condition. For 
Rd/Rb = 10, we always see the standard reentrant be- 
havior. For Z = 10, we can not see clear deviation from 
the VZK theory (compare curves a and b). With the de- 
crease of the magnitude of Z, the magnitude of Sas/a^ 
is enhanced much stronger (curves c and d). Although 
there are quantitative difference between four curves, the 
qualitative line shapes are similar to those predicted by 
the VZK theory. 



Rd/Rb=l 




FIG. 13: Sas/cTN is plotted as a function of T. Rd/Rb = 1 
and Erh = 0.01 Aq. a: VZK theory, b: Z = 10, c: Z = 1 and 
d: Z ^0. 

However, situation is different for decreasing Rd/Rb — 
1. For Rd/Rb ~ 1 the resulting Sas/a^ takes negative 
value and similar feature is also obtained for Z = 10 (see 



curves a and b). At the same time, for small magnitude 
of Z, 5as/aN takes positive value. As seen from these 
results, the deviation from VZK becomes significant for 
small magnitude of Z with Rd/Rb < 1. 



IV. CONCLUSIONS 

In the present paper, a detailed theoretical investiga- 
tion of the tunneling conductance of diffusive normal 
metal / conventional superconductor junctions is pre- 
sented. Even for conventional s- wave junctions, the inter- 
play between diffusive and interface scattering produces 
a wide variety of line shapes of the tunneling conduc- 
tance: ZBCP, ZBCD, U-shaped, and rounded bottom 
structures. There are several points which have been 
clarified in this paper. 

1. When the transparency of a junction is sufficiently 
low and Erh is small, the ZBCP appears for Rd < Rb- 
With an increase in Rd, the ZBCP changes into a zero 
bias conductance dip (ZBCD). For large Eth with the 
same order of Aq, the ZBCP is only expected for small 
values of Rd- With increasing Exh: the ZBCP vanishes 
when ETh > Aq is satisfied. For the low transparency 
limit, the results obtained by us are reduced to those 
by the VZK theory where the KL boundary condition is 
used. 

2. When the transparency of the junction is almost 
unity, aT{e.V) always have a ZBCD except for the case 
of vanishing Rd- 

3. The proximity effect can enhance (reduce) the tun- 
neling conductance of junctions with low (high) trans- 
parency. 

The above mentioned criteria for the existence of a 
ZBCP agree with available experimental data. However, 
we are not aware of an experimental observation of ZBCD 
in highly transparent junctions. 

In the present paper, the superconductor is restricted 
to be a conventional s-wave superconductor. However, 
it is well known that a ZBCP also appears in uncon- 
ventional superconductor junctions^S*^, the origin of 
which is the formation of niidgap Andreev bound states 
(MABS)^^. Indeed, a ZBCP has been reported in var- 
ious superconductors that have an anisotropic pairing 
symmetry4ii^ It should be remarked that the line shape 
of the ZBCP obtained in the present paper is quite differ- 
ent from that by MABS. In the present case, the height of 
aT{eV) never exceeds unity and its width is determined 
by ETh, while in the MABS case, the peak height is pro- 
portional to the inverse of the magnitude of Tm and the 
width is proportional to Tm- 

The proper theory of transport of unconventional junc- 
tions in the presence of MABS has been formulated^Si^ 
only under the conditions of ballistic transport. Recently, 
this theory has been revisited to account for diffusive 
transport in the normal metal in Refi^, where an ex- 
tension of the circuit theory was provided for unconven- 
tional superconductor junctions. A general relation was 
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derived for matrix current, B in Eq. (4), which is avail- 
able for unconventional superconductor junctions with 
MABS. An elaborated example demonstrated the inter- 
play of MABS and proximity effect in a d-wave junction. 
It is actually quite interesting to apply this novel cir- 
cuit theory to the actual calculation of arieV) as in the 
present paper. Such a direction of study is important 
in order to analyze recent tunneling experiments where 
mesoscopic interference effects were observed in high-Tc 
junction systems^. We will present the obtained results 
in the near future^. 

In the present study, we have focused on N/S junc- 
tions. The extension of Nazarov's circuit theory to long 
diffusive S/N/S junctions has been performed by Bezug- 
lyi et aP^. In S/N/S junctions, the mechanism of mul- 
tiple Andreev reflections produces the subharmonic gap 



structures on I-V curveo'^^i'^^i^^i'^^i^°i^^i'^^i^^ and the situ- 
ation becomes much more complex as compared to N/S 
junctions. Moreover, in S/N/S junctions with unconven- 
tional superconductors, MABS lead to the anomalous 
current-phase relation and temperature dependence of 
the Josephson current^. An interesting problem is an 
extension of the circuit theory to S/N/S junctions with 
unconventional superconductors. 

The authors appreciate useful and fruitful discussions 
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supported by the Core Research for Evolutional Science 
and Technology (CREST) of the Japan Science and Tech- 
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percomputer Center, Institute for Solid State Physics, 
University of Tokyo and the Computer Center. 
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